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Ph! ABSTRACT: Following recent developments, we employ the AdS/CFT correspondence 

to determine the drag force exerted on an external quark that moves through an M = 4 super- 
Yang-Mills plasma with a non-zero R-charge density (or, equivalently, a non-zero chemical 
potential). We find that the drag force is larger than in the case where the plasma is neutral, 

H ' but the dependence on the charge is non-monotonic. 
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1 Introduction and Summary 



Central collisions of gold nuclei at RHIC are believed to produce the long-sought quark 
gluon plasma (QGP). RHIC experiments have found evidence of strong collective behavior, 
and the consensus is that we are dealing with a strongly-interacting liquid that can be 
modelled by hydrodynamics (see, e.g., pQ for a recent review). The hydrodynamic regime 
is completely characterized by transport coefficients (shear and bulk viscosity, etc.). RHIC 
data suggests that the QGP viscosity should be fairly small — the ratio of shear viscosity to 
entropy density has been estimated [2] to be rj/s ~ 0.1 — 0.2 , although uncertainties in this 
value are still quite large. Unfortunately, a calculation of the hydrodynamic coefficients from 
first principles, i.e, from the underlying microscopic theory, is currently out of our reach in 
the strong-coupling regime. 

The gauge/gravity or AdS/CFT correspondence [3][I] has been used to investigate observ- 
ables in various interesting strongly-coupled gauge theories [5] where perturbation theory is 
not applicable. Policastro, Son and Starinets pioneered the study of hydrodynamic gauge 
theory properties using AdS/CFT [HIT]. In [6], these authors computed the shear viscosity 
f] of strongly-coupled Af=4 super- Yang-Mills (SYM) theory in 3 + 1 dimensions and at finite 
temperature. They found that the ratio of shear viscosity to entropy density equals 1/47T, a 
result that was later argued to be universal, in the sense that it applies to any gauge theory 
described by a supergravity dual, in the limit of large 't Hooft coupling [H]Jj These results 
raised the tantalizing possibility of using AdS/CFT to study the QGP. Transport coeffi- 
cients of different thermal gauge theories have been calculated in [7]- [22]. Work attempting 
to narrow the gap between the gauge/gravity duality and RHIC may be found in [231 El]- 

RHIC data also confirmed the existence of jet quenching in high-energy heavy ion collisions 
[25] and showed that the observed suppression of hadrons from fragmentation of hard partons 
is due to their interaction with the dense medium [26]. In a series of interesting and very 
recent papers, the AdS/CFT machinery has been brought to bear on the phenomenon of 
jet quenching and the associated parton energy loss. The authors of [27] suggested that 
the jet-quenching parameter q used as a measure of energy loss in phenomenological models 
[28] [29] could be defined non-perturbatively in terms of a lightlike Wilson loop, and then 
employed the gauge/gravity duality to compute this loop in strongly-coupled Af = 4 SYM. 
The calculation was extended to the non-conformal case in [30]. The works [31] [32] studied 
the dynamics of moving strings in the AdSs-SchwarzschildxS 5 background to determine the 
drag force that a hot neutral M = 4 SYM plasma exerts on a moving quark. In [33] , through 
an analysis of small string fluctuations, the authors determined the diffusion coefficient for 
the plasma, finding a result that is in agreement with [3T] . 

In this note we extend the computation of [31] [32] to the near-horizon geometry associ- 
ated with rotating near-extremal D3-branes, which allows us to examine how the drag force 

1 The leading order correction in inverse powers of the 't Hooft coupling was determined in [9]. 
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changes upon endowing the plasma with a non-zero charge (or, equivalently, chemical po- 
tential) under a U(l) subgroup of the £77(4) R-symmetry group. The fact that the QGP 
produced at RHIC is also found to carry a small but non-zero charge density (associated to 
baryon number) cannot but add interest to our calculation. But of course, since the string 
theory dual of QCD is not known, one cannot overemphasize that, for multiple reasons, com- 
parison of AdS/CFT results with real data is, in general, a risky business (for discussion, 
see [S1E2]). 

In Section [2] we review the relevant properties of the rotating black three-brane geometry, 
and discuss the embedding of the Nambu-Goto string in this background. Section [3] contains 
the calculation of the drag force. The dependence of our general result (I3.20p on the plasma 
temperature T and charge density J is shown in Fig. 1, obtained by solving the sixth- 
order equation (13.211) numerically. For small charge density, this equation can be solved 
perturbatively, and (13.241) gives the resulting drag force at next-to-leading order. Our main 
conclusion is that the force in the charged plasma case is larger than in the neutral case. 
Interestingly, we find the same result for the force in two cases where the external quark 
is coupled in different ways to the SYM scalar fields, despite the fact that in one ('polar') 
case the quark is neutral under the global U(1)r, while in the other ('equatorial') case it is 
charged. 

While the first version of this paper was in preparation an overlapping preprint [31] ap- 
peared. We will discuss its relation to our work at the end of Section [31 

2 String in a Rotating D3-brane Background 

The AdS / CFT correspondence [31 0] equates the physics of an external quark in a finite- 
temperature neutral M = 4 SYM plasma to that of a fundamental string that lives on the 
near-horizon geometry of a stack of static near-extremal D3-branes |35j . This connection 
was employed in [311 [32] to determine the drag force that the neutral plasma exerts on the 
quark. Our aim in this paper is to consider instead a plasma that is charged under the 
S77(4) R-symmetry, and so we must analyze a string that ploughs through a rotating D3- 
brane background. For simplicity we will restrict attention to the case where only one of the 
three S77(4) angular momenta is non-zero. The corresponding solution was obtained in [36] 
(see also [37J)p we will follow here the conventions of [39]. In the near-horizon limit, the 
metric is 



ds 2 



i 



{-hdt 2 + dx 2 ) + Vh 




(2.1) 



+r 2 (Ad6 2 + A sin 2 6d(j) 2 + cos 2 Od^l)] , 



2 The metric for the case with all three angular momenta turned on can be found in [38] . 
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where 

H 
A 



R 4 

r 4 A ' 

I 2 cos 2 e 



x I 2 
A = 1 + 



h = 1 



r 4 A ' 
1 / Z 2 r 4 



A y r 2 7 ' 
# 4 = 4vr^ s / 4 . (2.2) 

This geometry has an event horizon at the positive root of h(rn) = 0, 

rl = \ (\A 4 + 4r 4 - / 2 ) , (2.3) 

and an associated Hawking temperature, angular momentum density, and angular velocity 
at the horizon 

T =V^^<, J=7^Ws> Q = %> ( 2 - 4 ) 
2nR 2 rQ v 647r 4 g 2 Z° r£R 2 

which translate respectively into the temperature, R-charge density and R-charge chemical 
potential in the gauge theory. For J ^ 0, there is an ergosphere between r = ru and the 
stationary limit surface r = r s {9) defined by h(r s , 9) = 0. 

An external quark (a pointlike source in the fundamental representation of the SU(N) 
gauge group) that moves at constant velocity in the x l = x direction and carries a charge 
under the same global £7(1) r C SU (4) r as the plasma corresponds to a string whose embed- 
ding in the geometry (12. ip and in the static gauge a = r, r = t is of the general stationary 
form 

X\r,t) =vt + £(r), 9(r,t)=vt + C(r), 0(r, t) = ut + x(r), (2.5) 

with all other fields vanishing. As usual, the behavior at infinity describes the gauge theory 
in the extreme UV, where the pointlike quark is inserted, so the string must have a single 
endpoint on the boundary. The corresponding boundary conditions^ 

X(r,t) — >• Xoo(f) = vt, 9(r,t) — »• 9oo(t) = vt + Coo, 0(r,t) -> 0oo(t) = cut as r -> oo, 

(2.6) 

specify both the trajectory of the quark (which determines its coupling A At (A 00 (t))^A^ J (t) to 
the gluonic fields,) and its global U(1)r charge (which controls its coupling $ i (X 00 (t))nJ 50 (t) 



3 Notice that an appropriate choice of origin allows us to set £oo = = Xoc without loss of generality. 
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to the scalar fields, with n^t) G R 6 the unit vector which points towards the north pole 
of the S 3 C S 5 lying at polar angle 9oo(t) an d azimuthal angle 4>oo{t) [M])- The tail of the 
string is associated instead with a flux tube that could be mapped out with the methods of 
[42l |43]FI In the present case the tube which codifies a fixed gauge- and scalar- field pattern 
that follows the moving quark and becomes wider as one moves back along the — x direction 
at a given time (or, equivalently, as time elapses at a fixed location in x) [32] . 

The equations of motion for the string are obtained from the Nambu-Goto Lagrangian 



C = = -V- det(d a X^d b X-G^) , (2.7) 

where G^ u (/^, v — 0, . . . , 9) is the spacetime metric (12. ID and g ab the induced metric on the 
worldsheet. The stationary embedding (12.51) is such that d t (dC/dX^) = V/i, and as a 
result, the equations for the cyclic variables X(r, t) and </>(r, t) amount to the statement that 
the conjugate spatial momentum densities 

dC dC 

n * = n * = OX 1 ' ^ = ^ = dft ^ ' 

are constant. Our main goal is to determine the value of tt x , which controls the ^-component 
of the force that each segment of the string exerts on its larger-r neighbor 



F x = ir x . (2.9) 

In the situation of interest the string trails behind its boundary endpoint (which is being 
pulled in the +x direction by an external agent), so we should have F x < and therefore 
ir x < 0. 

The equation of motion for 6(r,t), on the other hand, is rather complicated, due to fact 
that £ has explicit ^-dependence. In particular, and in contrast with the non-rotating case, 
a constant-6 1 solution is possible for I ^ only if 9^ — 0, 7r or 9^ — ir/2, just like in the 
v = case studied previously in [ID] . For simplicity, we will restrict attention to these cases, 
which describe a string that lies respectively perpendicular and parallel to the rotation plane 
of the black brane. 

• In the 'polar' case 9^ = or tt, the string points at constant angles 8(r,t) = 0, it and 
0(r, t) = 0, and reaches the horizon at the pole, where the ergosphere has vanishing 
width (h = h, so the stationary limit surface coincides with the horizon). Starting 
from the definition of ir x in (1 2 . 8 j) . it is then easy to infer that 



X' 



G„ {-GujGxx - v 2 ) ^2 io) 

\ —Gu(—GttG xx — ix x ) 



4 After this work had appeared as a preprint on the arXiv, the relevant calculations were carried out both 
for the case of vanishing [44] and non-vanishing [45] R-charge density. 
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where, in accord with the previous discussion, the overall sign has been chosen in such 
a way that ir x < implies X' > 0, meaning that the string trails behind its boundary 
endpoint. 

In the 'equatorial' case 9^ = 7r/2, the string is at constant polar angle 9(r,t) = tt/2 but 
can have a nontrivial azimuthal profile </>(r, t) — out + x( r )- One can again invert (12.81) 
to obtain X' and 0' in terms of tt x , 7r^, v, to. To simplify the algebra, we will henceforth 
concentrate on the case u = 0, where the string does not rotate. This choice is clearly 
valid from the point of view of the gauge theory, where it amounts to specifying a 
time-independent coupling between the external quark and the scalar fields. On the 
other hand, it might seem counterintuitive from the gravity perspective, because we 
expect the string to penetrate beyond the ergosphere, where a point particle could not 
avoid rotating. Nevertheless, we will see in the next section that it is possible to find 
physical configurations where the string does precisely that. Setting u = 0, then, we 
find that the definitions ( 12. 81) imply 



X' 



D 



G r 



■K^GtyG^v + tt x {—G 2 t(j) + G^(G tt + G xx v 2 )} 



G rr G xx 



In^z-Gt^G^v + 7T^{— + G^(G tt + G xx v 2 )} 
+G xx {n 2 Gu — {G t( f, — G^Gtt){Gu + G xx v 2 )} , 



(2.11) 



where again the overall signs were chosen such that n x < for the desired trailing 
string. 



3 Drag Force in a Charged Plasma 

We are now ready to compute the drag force exerted by the charged plasma on the quark, 
following [2H The first step is to demand that the solutions to (12.101) and (12. lip be 
well-defined over the entire region rn < r < oo. This condition is nontrivial because the 
solutions involve square roots of quantities that can become negative. As we will now learn, 
insisting that this does not happen will fix the values of 7r x , 7r^. As in [3"T| I3"2], in the analysis 
to follow a special role will be played by the velocity-dependent radius r v {9) that is the 
largest root of the equation h(r, 6) — v 2 = 0, 



rv(9) 2 = - hpcos 4 # + r -^-Z 2 cos 2 # . (3.12) 
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In the polar case X' is given by equation (12.101) . The numerator involves the factor 
—G tt /G xx — v 2 = h — v 2 , which according to (13.121) changes sign at 



r,(0) 



\ 




(3.13) 



The only way we can prevent X' from becoming imaginary for r < r v (0) is by choosing 
a value of ir x such that the denominator also changes sign at r v (0), i.e. 



7T, 



h(r v ) 



R 2 VT 



(3.14) 



The string profile can then be determined by integrating the equation obtained by 
plugging (1533)1 into fl2~T0l) . 



vr\R 2 



r 4 + l 2 r 2 — r 4 



(3.15) 



In the equatorial case X' and 0' are given by (12. lip . Using the explicit metric com- 
ponents (12.11) . it is easy to convince oneself that, out of the various factors inside the 
square roots, only the function 



D(r) 



o 



+ 



i? 4 



r 2 + 



o 



(3.16) 



r 2 + 



I 2 77" 

; R 4 R 4 





'l-v 2 ' 


r 4 + 


R 4 



runs the risk of becoming negative. This function clearly approaches +oo both at 
r — > and r — * oo, so it must have at least one minimum at some intermediate point. 
Since r 2 D(r) is quartic in the variable r 2 , it could in general have as many as four 
distinct roots at positive values of r. Numerical calculation shows that, for generic 
values of 7r x ,7r^, the function D(r) has exactly two positive roots, r±,r2 > Th, and a 
single minimum in between. This implies that D(r) < (and the solution (12. lip is 
ill-defined) in the intermediate range r± < r < r2, an undesired feature that can only 
be avoided by choosing ir x , in such a way that the D(r) curve has its minimum at 
D(ri — r^) — 0. Analytically, this corresponds to demanding that the r° and r 4 terms 
in ( 13.16P vanish, implying that 



7I> 
7T* 



iVT 



R 2 VT^ 



(3.17) 



6 



after which D is found to take the form 



D(r) = -g;- 



4 

r 4 r 



V 2 



(3.18) 



This function is manifestly non-negative, and evidently has a minimum and a double 
root at r /(l — f 2 ) 1 / 4 . Notice that this is none other than the critical radius r v (ir/2) 
defined in f l3TT2|) . 

The string profile in this case follows from integration of the equations of motion (12. lip , 
which simplify drastically after use of (13.171) : 



r 4 + i z r z — Tq 
lr 2 



r 4 + l 2 r 2 — r 4 



These expressions are clearly well-defined over the entire range m < r < oo. One 
might have also worried that the string worldsheet could become spacelike at some 
point inside the ergosphere, but in fact (I3.19P leads to y/—g = yl — v 2 , which is 
manifestly real. So as promised, we have been able to find a physical configuration 
where the string penetrates into the ergosphere and yet does not rotated What has 
happened is that the string has managed to adopt an azimuthal profile <p(r) that leads 
to a precise balance between the string tensile force and the inertial drag due to the 
rotating geometry. Notice that, according to (13.171) . this requires an external agent to 
pull on the boundary endpoint of the string, exerting the non-zero force n^. Notice 
also that the sign of 0' in (13.191) correctly reflects the fact that, as r increases, the 
string should wind in the direction opposite to the sense of rotation. 

The second step is to use (13. 14|) or (I3.17P in (12. 9p . to read off the force with which the 
string pulls back on its boundary endpoint, which by the AdS/CFT correspondence should 
be identified with the drag force that the charged plasma exerts on the moving quark. 
Surprisingly, the result is the same in the polar and equatorial cases, 

„ dp x r 2 /R 2 v , . 

F -i = -W7r=7' (3 ' 20) 

We learn then that the value of the U(1)r charge of the external quark affects the gluonic 
and scalar field distributions set up by the quark (which, as explained in Section [21 are 
encoded in the shape of the string tail (13.151) or (13.191) ). but does not alter the total drag 



5 In the process, we have also learned that the solution ran the risk of being ill-defined not just inside 
the ergosphere, but in the region below the critical radius r v {9) defined in (|3.12|) . which coincides with the 
stationary limit surface r s (6) only for v = 0. 
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force experienced by the quark. Notice that the end result (13.201) is identical to the one 
given in equation (12) of [32], except for the fact that r# has been replaced here by r®. The 
two radii indeed agree for the non-rotating case, but are related through (I2.3P in the general 
case, which shows that the drag force does have an interesting dependence on the R-charge 
density of the plasma. 

The third and final step is to rewrite (13.201) in terms of gauge theory parameters. Using 
(12. 3D and (12. 4ft one finds that p = r^/W^I^T 4 satisfies the sixth-order equation 



lQp b - p J 



Ac 2 p 4 + 8c A p 3 -c 6 p + c 8 







where we have defined a dimensionless charge parameter 

J 



2nN 2 T 3 



(3.21) 



(3.22) 



In the neutral case c = 0, equation (13.211) implies that p = 1/16, from which one immedi- 
ately recovers the well-known relation r = ttR 2 T. For small charge densities one can obtain 
ro by solving (I3.21j) in an expansion in powers of c. The result at next-to-leading order is 



r 



ttR 2 T 



1 + 



4J 2 



8rrl2\ 



7T 2 iV 4 T 6 



0(jyN s T 



(3.23) 



which together with (12.21) . (I3.20p . p x = mv/^/1 — v 2 and g\ M = ^ixg s leads to the final result 

8J 2 



dPx 
dt 



71 



1 + 



7T 



2jv 4 T 6 



(3.24) 



From this we can read off the exponential relaxation time 



To 



Ji2 



8J 2 



7r 2 iV 4 T 6 



0(J*/N*T 



(3.25) 



We learn here that turning on a nonzero R-charge density for the plasma increases the drag 
force exerted on the heavy quark (or equivalently, decreases its relaxation time). Notice that 
the J-dependent terms involve inverse powers of N, as has been found when computing other 
properties of the charged plasma (e.g., its entropy density [39]). 

The behavior of the drag force for larger charge densities can be determined by solving 
(I3.2ip numerically. The result is shown in Fig. [Tjj which displays the ratio between the force 
in the charged case and that in the neutral case. Real solutions exist only in the finite range 
< c < c max ~ 0.2349. The main feature is the maximum at Cq ~ 0.199, where the force 
ratio reaches a value ~ 1.299, and beyond which it decreases to ~ 1.184 at c max . 



6 The original preprint that we posted to the arXiv gave the arbitrary-charge result only implicitly, through 
(|3.21|) . We have added the plot to this later version in order to facilitate comparison with [3U [57] . 
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Figure 1: Drag force as a function of the dimensionless charge parameter c = J/2irN 2 T 3 . 

As explained in the Introduction, the authors of [27] advocated the use of a lightlike 
Wilson loop to define a parameter q meant to be an alternative characterization of energy 
loss in a thermal plasma. The works [581 ESI E7J [60] followed this procedure to determine q 
for the charged M = 4 SYM plasma, obtaining results that are rather similar to the drag 
force computed above. The drag force curve given in Fig. [1] is in fact nearly identical to 
the analogous q curve depicted (as a function of £ = 4y / 2c) in Fig. 2 of [57] . (Particularly 
striking is the fact that the two ratios coincide at £ max = A\/2c max .) This suggests that, at 
least in the current setting, both quantities indeed capture essentially the same physics. 

It is also interesting to compare our results with those obtained in the closely related 
work [31], which appeared while the first version of this paper was in preparation. After 
carrying out a general drag force calculation applicable to any gauge theory whose dual 
gravity description involves an asymptotically AdS^+i geometry, the author of [31] considers 
as a specific example the case of d = 4 R-charged M = 4 SYM plasma, with all three of the 
independent charge parameters (corresponding to the independent rotations within SU(4)r) 
turned on. In the present paper we have examined precisely this system, in the special 
subcase where two of these parameters are taken to vanish. 

There is, however, an important difference in our calculations: whereas the string stud- 
ied here moves in the ten-dimensional background ( 12. ip . which is a solution of Type IIB 
supergravity describing the near-horizon geometry of a stack of rotating D3-branes, the 
string analyzed in [34] lives on the five- dimensional 'STU' background, which corresponds 
to a charged black hole solution of jV = 8 D = 5 gauged supergravity [EE]. The latter 
theory is believed to be a consistent truncation of the Kaluza-Klein S 5 reduction of Type 
IIB supergravity, and in [37] it was shown that the STU and spinning D3 backgrounds are 
indeed related in this manner: the background described by equations (5.1)-(5.5) of [31] is 
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shown there to be a non-trivial truncation of our (12.11) . with the identifications n% = {I/th) 2 , 
&2 = ^3 = 0. The drag force deduced from the STU background in this case is given in equa- 
tion (5.12) of [3l], which is then to be compared with our spinning D3 result (I3.20p . The 
two forces are clearly very different: the result of |34j displays a complicated dependence on 
the velocity v and the charge parameter k = (l/rn) 2 , which is completely unlike the simple 
dependence seen in (13.201) . 

This discrepancy seems to suggest that: i) the putative uplift from five to ten dimensions 
of the string configuration considered in |34j would have a complicated (but presumably 
still stationary, as in ( 12.51) .) radial dependence for both the polar and the azimuthal angles, 
unlike the polar and equatorial profiles considered here; ii) string configurations more general 
than the ones considered here (e.g., a non-polar string dual to a quark whose coupling 
with the SYM scalar fields is time-independent) could potentially experience a drag force 
whose functional dependence on the relevant parameters is much more complicated than 
the one seen in (I3.20p . It would be interesting to explore this second point by searching 
for explicit configurations of this type in ten dimensions. The first point is obscured by the 
non-trivial nature of the truncation that connects the two backgrounds. In particular, given 
that the string of [3l] does not couple to the five- dimensional gauge field A t & = —lrl/R 2 r 2 A, 
which descends from the ten-dimensional metric component Gt * [37J , it would seem like the 
angular dependence of its uplift should be special enough to somehow decouple the angular 
problem from the five- dimensional geometry, which is all that is captured by the Nambu- 
Goto action used in [HI]. Our polar string seems to do precisely that, so it is not clear why 
this configuration appears not to be accessible from the five-dimensional perspective. 

A second possibility is that the string of [HI] is in some sense smeared over the S 5 . This 
would imply that its dual quark is qualitatively different from the one considered here, 
because its coupling to the SYM scalar fields would be somehow averaged over. A third 
possibility that should be kept in mind is that, even though the categorization of Af = 8 
D = 5 supergravity as a consistent truncation of Type IIB supergravity would guarantee 
that any solution of the former can be uplifted to a solution of the latter, it could turn out 
not to be possible to capture the dynamics of a string that lives in the uplifted background 
in terms of a Nambu-Goto action that senses only the five-dimensional geometry. 

The results obtained in [HH [321 [HI] and the present paper refer to a colored object, the 
external quark. For comparison, it might be worth computing the drag force on a color- 
neutral object, such as a mesorjll [l£l HH [35] or a baryon [271 HH1 EH]- We hope to return to 
this and related problems in future work. 

7 After this work had appeared as a preprint on the arXiv, the corresponding calculations were carried 
out in [5(2 E21 EI] , and related work was reported in ESI EH EH [56] . 
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